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> \ \ In the paper "Big Bang Nucleosynthesis and Active-Sterile Neutrino Mix- 

O \ ing: Evidence for Maximal Muon-Neutrino/Sterile-Neutrino Mixing in Super 



W | Kamiokande" (http:/ /xxx. lanl.gov/abs /astro-ph/9810075| ) , we suggested that 



to evade the Big Bang Nucleosynthesis exclusion of the muon neutrino to 
sterile neutrino oscillation explanation of the Super Kamiokande data, the 
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^C) , tau neutrino must have a mass over about 15 eV and it must mix with a 

lighter sterile neutrino. A stable tau neutrino with this mass is inconsis- 

00 

\ tent with cosmological structure formation. In a comment on our paper 



( http: / / xxx.lanl.gov / abs / astro-ph /981 1067|) , Foot and Volkas argued that our 



43 

result is incorrect and that the required tau neutrino mass should be much 
lower. Here we back up our original result with a more detailed calculation. 
We show that the argument of Foot and Volkas is invalid, most likely due 
to an insufficient energy resolution in the low energy part of the neutrino 
spectrum. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35. +c; 95.30.Cq 
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The Super Kamiokande atmospheric neutrino data have provided evidence for muon 
neutrino oscillation as a result of a maximal or near maximal mixing between and either 
v T or a sterile neutrino v s The mass-squared- difference involved in the mixing is ~ 10 -3 
eV 2 The maximal or near maximal z/ M <-> v s mixing explanation, however, would violate 
the Big Bang Nucleosynthesis (BBN) bound by nearly completely equilibrating a fourth 
neutrino flavor v s in the early universe and, hence, increasing the primordial 4 He yield 
above the observed level One mechanism proposed to evade this bound is to invoke 

a large lepton number asymmetry (> 1CT 5 ) at the BBN epoch which acts to suppress the 

<-» v s transformation by matter effects. A clever scheme along these lines suggested by 
Foot and Volkas |4j] is to have a massive v T mix with a lighter sterile neutrino v s i (which 
could in principle be the same v s ). The resonant transformation of tau neutrinos to sterile 
neutrinos via matter-enhanced mixing at the BBN epoch would generate a tau lepton number 
asymmetry L Vt that grows with time @||. The v T «-> v s < resonant transformation must 
occur before any significant <-> v s transformation can occur, so that the L Vt generated 
can subsequently suppress the <-> v s transformation. However, one twist in this scheme 
is that while L Vt grows, it crosses a parameter region where the MSW (Mikheyev-Smirnov- 
Wolfenstein) matter effect causes either (if L Vr > 0) or (if L Ut < 0) to resonantly 
transform into v s or its antiparticle The result is a newly generated L Vfil with a 

sign opposite to that of L Vt . This L„ acts to counter the suppression effect of L Vt on the 
Vp, *-* v s transformation. This is because the matter-antimatter asymmetry contribution to 
the effective potential of the <-> v a system is oc 2L V + L Vr [0. The countering effect 
of L v is immaterial only if L Vt is sufficiently large at the time of the or resonant 
conversion. In turn, this requires the mass-squared-difference of the v T <-> v s > mixing to 
satisfy m 2 T — m 2 ; > 300 eV 2 , which implies m UT > 15 eV [0. Tau neutrinos this massive 
are incompatible with cosmological structure formation, and therefore would have to be 
unstable 0. 

In papers and in a comment to our aforementioned paper, Foot et aJ. investigated a 
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similar problem and argued for a much lower limit m 2 l>r — m? u / > 16 eV 2 0. The difference 
between the required threshold v T mass in our calculations and theirs may stem from differ- 
ences in the energy resolution employed. Namely, we believe that an accurate treatment of 
the counter-suppression effect of L v requires a very fine energy resolution in the low energy 
part of the mu neutrino spectrum. 

We agree with Foot and Volkas that L Vt grows exponentially at the initial stage of 
lepton number generation but then subsequently approaches an asymptotic T" 4 growth 
(where T is the temperature of the universe). Foot and Volkas argued that since most of 
V/j, or underwent resonances during the exponential growth phase, there was not enough 
time for any significant generation of L v . Our results, on the other hand, show that the 
growth rate of L Vt is not the central issue in the problem. Rather, since L Ut itself is small, 
especially during its initial phase of exponential growth (e.g., < 10~ 8 ), even the resonant 
conversion of a tiny fraction of or into sterile neutrinos may generate a competing and 
significant L v « —L Ut /2. Once such an L Vfi has been generated, the potential responsible 
for the <-> v s transformation is driven close to zero, rendering the suppression from 
L Vt ineffective. Therefore, in order to suppress the <-> v s mixing required for Super 
Kamiokande, the v T «-> v s i mixing must generate an L Vt that is much larger than its induced 
L v at any moment. 

There are two resonances involved in the problem: the v T <-> v a > resonance that generates 
L Ut ; and the L,, r -induced <-> u s resonance that generates a competing L u ^. Because 
the effective potentials of active-sterile neutrino mixings are neutrino energy dependent, at 
any given temperature each resonance occurs only in a narrow energy bin in the neutrino 
energy spectra. As the temperature of the universe drops, a particular resonance energy bin 
gradually moves to higher neutrino energy, eventually sweeping across the entire neutrino 
energy spectrum. The narrow widths of the resonance energy bins and the rapid decrease of 
the effective mixing angles outside the resonance energy bins enable a simple analytical/ semi- 
analytical calculation. 

In our semi-analytical numerical approach, we track the neutrino transformations in 
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the resonant parts of the neutrino energy spectra. Again, these resonant regions in energy 
space consist of only narrow energy bins and are functions of the temperature and the 
lepton number asymmetry. This semi-analytical approach offers two distinctive advantages: 
a very high energy resolution of the neutrino spectrum (e.g.,AE/T < 10~ 6 ); and ease 
in understanding the physical processes involved, especially the interplay between the two 
coupled mixing systems. A high energy resolution is essential because the initial exponential 
growth of L Vt comes from minuscule differences between the resonance energies of the v T <-> 
v s i system and the v T «-> v s > system, and the tiny L Vt at this stage can be easily matched 
by a competing L v . Energy bins that are too coarse may therefore result in an incorrect 
account of the L Ut growth and an omission of the counter-balancing effect of L v when L Ut 
is small. 

Figure 1 shows our result based on the same mixing parameters assumed in the figure of 
Foot and Volkas' comment ( |astro-ph / 98 1 10671 ) to our paper. For the purpose of comparison, 
no simultaneous <-> v s transformation is assumed in this figure. Figure 1 is for the 
most part similar to that of Foot and Volkas. There are minor differences that are readily 
identifiable: (1) our result tracks T -4 more closely in the "power-law growth" epoch; (2) L Ut 
in our results does not switch sign at the initial point of growth. The sign difference is not 
surprising because of the chaotic character of the growth, which introduces a sign ambiguity 
to the problem ||. 

As an illustration, also plotted in Figure 1, are the \L Vt \ required for the <-> v s or 
<-> v s resonance to occur at or energies e res = p^l/T = 0.01, 1, 10. (The parameters 
for the i/fj, <-* v s mixing are Sm 2 = 10~ 3 eV 2 and sin26 l = 1.) It can be seen from the figure 
that the lower energy component of the or neutrinos encounters the resonance first 
when \L Vt \ is very small, and the resonance region moves through the or spectrum to 
higher neutrino energies as \L Vr \ becomes larger. 

For the potential proportional to L Ut + 2L U ^ to successfully suppress the *r+ v s trans- 
formation, the L Vfi generated by the or resonance has to be much smaller than L Ut in 
magnitude at any temperature as the or resonance sweeps through the entire spectrum, 
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f(e TCS )5e res R ^ < \ (l Vt + 2L Vfl ) (1) 

^res "J 

for any e res . In the equation, / is the Fermi-Dirac distribution function. 5e res is the energy 
width of the resonance. f(e Te8 )8e ies is therefore the fraction of mu neutrinos in resonance. R 
is the resonant transition rate. |5e res /e res | is the duration of the resonance at e res . The energy 
width of the resonance depends on whether the resonant transition is collision-dominated 
(with the quantum damping rate D ~ 0.5GfiT 5 e res > V x ) or oscillation-dominated (D < V x ): 

\ 2\Dde rcs /dV z \ i£D>V x 
Se ICS ~ . (2) 

[2\V x de Ies /dV z \ ifD<V x 
So does the resonant transition rate: 

( V x 2 /D if D > V x 
R^l . (3) 

\V X if D < V x 
The effective potentials for the <->• v s transformation are 

12 2 1 

V x = " J sin 20, V z « 22G|T 5 e±0.35G F T 3 (L i , T + 2L U/i ), (4) 

where GV is the Fermi constant, the "+" ("— ") sign is for z/ M (u^), and we employ natural 
units. (Here V y = 0.) The or resonance occurs at an energy e rcs ~ (L Ut +2L Uii ) /63G fT 2 , 
with de ves /dV z » (22GIT 5 )" 1 . The temperature T rcs at which the resonance occurs is almost 
independent of e res : T res « 22[(m^ T — m^J/leV 2 ]. (T res is an insensitive function of the 
v T «-> z/ s / vacuum mixing angle.) Furthermore, we only consider the case where e res <C 1. 
This is when L Ut is in its initial stage of growth (L Vt <C 10~ 7 ) and is most easily matched 
by a competing L„ . Given the small e res we have /(e res ) ~ e;? cs /1.8. Eq. (P can then be 
rewritten in the form 



12 2 12 

\m„ — mt 

1600G^T/ Cs |e res | 



< 80G F T r 2 es . (5) 



if D > V x , or 
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lOUUUfl res t res I e res | 



< 80GVT 2 . 



(6) 



if D < V x . 

We can further rewrite |e res | = -ffe r es|dlne re s/dlnT| where H = — dlnT/dt ~ 5.5T 2 /m p / 
is the Hubble expansion rate. The Planck mass is m p i « 1.22 x 10 28 eV. Then Eq. ([S]) 
becomes 



'ml T -m 2 



11/6 



leV^ 



> 2 x 10 4 e r -i 



din e r 



dlnT 



mr 



mr 



10- 3 eV 2 



(7) 



for D > V x . And Eq. 
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becomes 

2 \ !7/6 

-mi \ 



leV 2 



> 10 3 e" 3 

^ " '-res 



din e r 



dlnT 



2 2 

mf, — mi 



10" 3 eV 2 



(8) 



for D < V x . 

The value of |dlne res /dlnT| is related to the growth of L Ut by |dlne res /dlnT| « 
| d In L Ut /dlnT — 2| ~ |d In L UT /d lnT| (with safely ignored). Figure 2 shows 
| d In L Ut /dlnT | as a function of the v T <-> v s > vacuum mixing parameters, m 2 T — m 2 s and 
sin 2 26 1 ', in the initial exponential stage of L Vt growth when L Vr is <C 10~ 7 . |dlnLj, T /dlnT| 
can be approximately fit as 

d In L Vt 



dlnT 



6 x 10 6 sin2fl', 



(9) 



and is insensitive to m 2 ^ — m 2 s . The vacuum mixing parameter sin 29' of the v T <-> zv mixing 
must satisfy the BBN bound (m 2 ^ — m 2 J sin 4 20' < 10~ 7 eV 2 (the cut-off at the upper-right 
corner of Figure 2) 0. 

Eq. (0) , (pj) and (|]) show that the most stringent requirement on m 2 ^ — m 2 ; does indeed 
come not from v s resonances at e res ~ 3, but from resonances centered at the smallest 
possible e res as long as the or z/ M transition probability in that resonance energy bin is 
<C 1. This condition, expressed as 



R 



Se r 



< 0.1, 



(10) 
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can be rewritten as 



/ 2 2 \ -1/2 

/ m 2 VT - mi sl \ 


d In e res 


-1/3 


2 2 

mf, — rat 


2/3 


\ leV 2 ) 


dlnT 




10- 3 eV 2 





(11) 



regardless of the value of D/V x . It can be shown that e res is in the oscillation-dominated 
regime if 



e rP , < 0.25 



m,. 



ra; 



10-3 eV z 



1/2/2 2 

\ leV 2 



-1/2 



(12) 



mf, . comes from the oscillation- 



Therefore, the most stringent requirement on ra 2 , 
dominated regime for sin 2 29' > 10 -8 (while |dlne res /dlnT| > 10 3 ), and from collision- 
dominated regime for sin 2 29' < 10~ 8 (while |dlne res /dlnT| < 10 3 ). 

Combining Eq. (0) or (U) , Eq. (|) and Eq. ([TTD yields a requirement on the mass-squared- 
difference necessary to effect suppression of «-» z/ s transformation at the Super Kamiokande 
level: 

r2 Y /4 eV 2 for sin 2 20' > 10" 8 

(13) 



m VT - ra > 



200 (|m 2 M - m 2 J/10- 3 eV' 



(sin 20 



A-l/2 



2 , - ra 2 |/10- 3 eV 2 ) eV 2 for sin 2 29' < 10" 



This is in agreement with our previous work. 

Perhaps because of insufficient energy resolution at low energies, Foot and Volkas may 
have missed the crucial effect of small e res . Neutrinos with energies ~ 1% of the temperature 
are an insignificant fraction (~ 10~ 6 ) of the overall neutrino number. However, in this 
problem, they are the driving force which frees the «-> v s transformation process. This 
is simply because the suppressing lepton asymmetry from u T in this case is itself minuscule 
(e.g.,< 10- 7 ). 
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Figure Captions: 

Figure 1. The growth of the tau neutrino asymmetry as a result of the tau neutrino-sterile 
neutrino mixing, assuming m 2 ^ — m 2 / = 50 eV 2 and sin 2 29' = 10~ 4 . The intersections 
between the growth curve of L Vt and the dashed lines indicate when resonances occur for 
(if L Vt > 0) or (if L Ut < 0) neutrinos with momentum p. 

Figure 2. The initial rate of L Vt growth, dlnL/dhiT, as a function of the vacuum tau 
neutrino-sterile neutrino mixing parameters. 
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